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It  is  shown  that  in  a  spin-torque  microwave  oscillator  based  on  a  magnetic  nanocontact,  the  nature  of  the 
microwave  spin  wave  mode  generated  at  the  threshold  critically  depends  on  the  angle  between  the  external  bias 
magnetic  field  and  the  plane  of  the  free  layer.  When  the  external  bias  held  is  rotating  from  normal  to  in-plane 
orientation,  an  abrupt  transition  from  a  propagating  cylindrical  wave  with  the  frequency  higher  than  the 
frequency  of  the  linear  ferromagnetic  resonance  (FMR)  to  a  self-localized  standing  nonlinear  spin  wave 
“bullet”  with  the  frequency  lower  than  the  FMR  frequency  takes  place  at  a  certain  intermediate  angle  6cr  This 
transition  manifests  itself  as  an  abrupt  jump  (of  the  order  of  several  gigahertz)  in  the  generated  microwave 
frequency.  This  mechanism  of  mode  switching  might  explain  abrupt  jumps  of  the  generated  microwave  fre¬ 
quency  observed  in  recent  experiments  on  spin-torque  oscillators. 
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I.  INTRODUCTION 

It  was  predicted  theoretically1'2  and  confirmed  experi¬ 
mentally3-8  that  spin-polarized  current  passing  through  a  thin 
magnetic  layer  (“free”  layer  of  a  magnetic  layered  structure) 
can  excite  microwave  magnetization  oscillations  in  this 
layer.  Many  experimentally  observed  features  of  the  phe¬ 
nomenon  were  explained  in  a  series  of  theoretical  papers.9-12 
Although  it  is  clear  from  experiments  that  the  excitation  of 
microwave  oscillations  has  a  threshold  character  (i.e.,  it  is 
observed  for  sufficiently  large  currents  />  /[h)  and  that  the 
excitation  frequency  is  close  to  the  natural  ferromagnetic 
resonance  (FMR)  frequency  of  the  system,  the  exact  nature 
of  the  dynamic  spin  wave  modes  excited  at  the  threshold  was 
not  determined.  This  is  especially  true  for  the  case  of  cur- 
rent-driven  nanocontacts,3-6  where  magnetic  layers  are  not 
bound  in  plane,  and  there  are  no  reflective  lateral  boundaries 
that  could  create  an  effective  magnetic  resonator  (like  in  the 
case  of  nanopillars  with  well-defined  lateral  sizes7). 

The  dependence  of  the  microwave  excitation  in  a  mag¬ 
netic  nanopillar  on  the  angle  between  the  in-plane  magneti¬ 
zation  of  the  fixed  (reference)  layer  and  the  in-plane  direc¬ 
tion  of  the  external  bias  magnetic  field  in  the  free  layer  was 
studied  numerically  in  Ref.  13.  In  that  paper,  the  variation  of 
the  in-plane  orientation  of  the  bias  field  did  not  lead  to  any 
qualitative  changes  in  the  excited  spin  wave  modes,  but  sig¬ 
nificantly  influenced  the  transient  temporal  interval  during 
which  the  system  reached  the  steady-state  microwave  gen¬ 
eration. 

In  contrast,  in  the  case  of  the  nanocontact  geometry  and 
out-of-plane  orientation  of  the  external  bias  magnetic  field, 
the  variations  in  the  magnetization  angle  can  lead  to  a  quali¬ 
tative  change  in  the  nature  of  the  excited  spin  wave  modes. 

The  theoretical  analysis  of  spin  wave  excitation  by  spin- 
polarized  current  in  a  nanocontact  geometry  has  been  per¬ 
formed  in  two  limiting  cases:  in  the  case  of  normally914  and 


tangentially  (in  plane15)  magnetized  free  magnetic  layer.  It 
turns  out  that  the  nature  of  the  excited  spin  wave  modes  in 
these  two  cases  is  qualitatively  different. 

In  the  case  of  a  normally  magnetized  film,  both  linear9 
and  nonlinear14  analyses  predict  that  the  spin  wave  mode 
excited  by  spin-polarized  current  is  a  propagating  cylindrical 
spin  wave  with  the  wave  vector  kc—  1 .2/ Rc,  where  Rc  is  the 
nanocontact  radius.  The  coefficient  of  a  nonlinear  frequency 
shift  N  in  this  geometry  is  positive  N  >  0  (see  Ref.  11  for 
details),  so  that  with  the  increase  of  the  bias  current  I  (and 
the  amplitude  a  of  the  excited  spin  wave  mode),  the  fre¬ 
quency  of  the  excited  spin  wave  mode  also  increases.  Thus, 
the  frequency  a>(k )  of  this  mode  is  always  larger  than  the 
frequency  oj0  of  the  FMR  in  the  free  magnetic  layer: 

a>(k)  =  a>0  +  Dk2  +  \N\\a\2 ,  (1) 

where  k  is  the  wave  number  of  the  excited  mode,  oj0  is  the 
FMR  frequency,  and  D  is  the  spin  wave  dispersion  coeffi¬ 
cient  determined  by  the  exchange  interaction.  The  main  con¬ 
tribution  into  the  threshold  current  /[“  of  this  propagating 
mode  [see  Ref.  9  and  Eq.  (2)  in  Ref.  15]  is  provided  by  the 
radiation  losses  proportional  to  the  dispersion  coefficient  D. 
Therefore,  the  threshold  current  /[["  in  the  case  of  normal 
magnetization  does  not  vanish  even  in  the  case  of  an  arbi¬ 
trarily  small  damping. 

In  the  opposite  limiting  case  of  an  in-plane  magnetized 
nanocontact,  theoretical  analysis  performed  in  Ref.  15  dem¬ 
onstrated  an  absolutely  different  picture.  In  this  case,  the 
coefficient  of  the  nonlinear  frequency  shift  N  is  negative,  N 
<  0,  and,  therefore,  has  the  sign  that  is  opposite  to  the  sign  of 
the  spin  wave  dispersion  coefficient  D  (ND  <  0).  Thus,  apart 
from  a  linear  propagating  spin  wave,  the  in-plane  magnetized 
free  magnetic  layer  can  support  a  strongly  localized  non¬ 
propagating  spin  wave  mode  of  a  solitonic  type — spin  wave 
“bullet.”15  Although  it  is  well  known  that  two-dimensional 
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solitonic  wave  packets  (bullets)  in  conservative  systems  are, 
generally,  unstable  (see,  e.g.,  Refs.  16  and  17),  the  interplay 
between  the  spin-polarized  current,  creating  negative  damp¬ 
ing,  and  natural  positive  magnetic  damping  of  the  free  layer 
can  stabilize  a  spin  wave  bullet  of  a  certain  amplitude  a0  and 
size  €. 15,18  The  frequency  oj  and  the  spatial  size  €  of  the 
stable  spin  wave  bullet  excited  by  spin-polarized  current  are 
given  by  Eq.  (6)  from  Ref.  15, 


a)=(o0-\N\\a0\2,  €  =  ^  ^ ,  (2) 

Kl 

where  w0  is  again  the  FMR  frequency  in  the  free  layer. 

Since  the  frequency  of  the  spin  wave  bullet  (2)  is  shifted 
by  nonlinearity  below  the  spectrum  of  propagating  linear 
spin  waves,  the  bullet  mode  has  an  evanescent  nonpropagat¬ 
ing  character,  and  the  radiation  losses  for  the  self-localized 
bullet  vanish.  As  a  result  of  this  self-localization,  the  nonlin¬ 
ear  spin  wave  bullet  has  a  significantly  lower  excitation 
threshold  /^hul  than  a  linear  propagating  mode  [see  Eq.  (13)  in 
Ref.  15  for  details]. 

In  this  paper,  we  investigate  theoretically  the  nature  of 
spin  wave  modes  excited  by  spin-polarized  current  in  a  mag¬ 
netic  nanocontact  magnetized  at  an  arbitrary  out-of-plane 
angle.  We  assume  that  the  angle  0ml  between  the  direction  of 
static  magnetization  and  the  plane  of  the  free  magnetic  layer 
in  a  nanocontact  is  varied  continuously  from  0  (correspond¬ 
ing  to  the  case  of  in-plane  magnetization)  to  90°. 

We  demonstrate  that  the  nonlinear  nonpropagating  bullet 
mode  exists  in  a  wide  range  of  magnetization  angles  0 
=£  0int«  #lin>  where  6?lin  is  the  “linear”  magnetization  angle  at 
which  the  nonlinear  frequency  shift  vanishes,  N=  0.  We  also 
show  that  below  a  certain  critical  angle  0cr,  the  threshold  bias 
current  7th,  necessary  for  the  excitation  of  the  bullet  mode,  is 
lower  than  the  threshold  of  excitation  of  the  linear  propagat¬ 
ing  spin  wave  mode.  In  this  range  of  magnetization  angles 
0^  flintier,  the  spin  wave  mode  excited  at  the  threshold 
will  be  the  bullet  mode,  with  the  frequency  lower  than  the 
FMR  frequency,  while  for  larger  magnetization  angles  0CI 
=£  0ml=s  -77-/2,  the  spin  wave  mode  excited  at  the  threshold  is 
a  linear  propagating  spin  wave  mode.  The  dependence  of  the 
excitation  threshold  /th  on  the  magnetization  angle  0int  is  con¬ 
tinuous,  but  demonstrates  a  kink  at  the  transition  between  the 
two  excited  modes.  In  contrast,  the  angular  dependence  of 
the  frequency  oj(h  generated  at  the  threshold  has  an  abrupt 
jump  of  the  order  of  several  gigahertz  at  0int=  0CT  in  the  tran¬ 
sition  region. 


II.  ANALYTICAL  MODEL  OF  CURRENT-INDUCED 
MICROWAVE  GENERATION 


\a\2  =  (M0-M:)/2M0,  (3) 


where  M0  is  the  length  of  the  magnetization  vector  in  the 
free  layer  and  M;  is  the  projection  of  this  vector  on  the  axis 
z  defining  the  direction  of  static  magnetization  in  the  free 
layer.  This  equation  can  be  written  in  the  form15 


da 

dt 


■  i(u> o  -  DA  +  N\a\2)a  -  To  +  cr//|  —  j(l  -  |a|2)a, 


(4) 


where  oj0  is  the  frequency  of  linear  ferromagnetic  resonance, 

<W0  =  \  Bfl(%  +  U>M  cos“  $int)  .  (5) 

yHmv  ojm  =  AttjM {],  y~2.8  MHz/Oe  is  the  gyromag- 
netic  ratio,  D  =  (oM\lx(d(i>0/ dwH)  is  the  dispersion  coefficient 
of  spin  waves  due  to  inhomogeneous  exchange  interaction 
(Xex~5  nm  is  the  exchange  length),  A  is  the  two-dimen¬ 
sional  Laplace  operator  in  the  film  plane,  N  is  the  nonlinear 
frequency  coefficient  given  by 

l),  (6) 

E  =  aG(x)0(d(oQ/ f)(oH)  is  the  linear  Gilbert  damping  rate,  cr 
=  eg/bLB/2eM0dS  (e  is  the  spin-polarization  efficiency  de¬ 
fined  in  Refs.  1  and  9,  g  is  the  spectroscopic  Lande  factor,  /jlb 
is  the  Bohr  magneton,  e  is  the  modulus  of  electron  charge,  d 
is  the  thickness  of  the  film,  and  S=ttR2  is  the  area  of  the 
nanocontact),  I  is  the  charge  current  through  the  nanocontact, 
and  f(x )  is  the  dimensionless  function  describing  current  dis¬ 
tribution  across  the  nanocontact.  In  this  work,  we  will  as¬ 
sume  uniform  current  distribution,  i.e.,  fix )  =  I  for  x<  1  and 
f{x)  =  0  otherwise. 

In  Eqs.  (4)-(6),  Hm  and  9ml  are  the  interna!  bias  magnetic 
field  magnitude  and  out-of-plane  angle,  respectively,  con¬ 
nected  with  the  external  values  Hexl  and  f?ext  by  the  usual 
electrodynamic  boundary  conditions: 


^ext  ^OS  0ext  —  Hmt  COS  0mt, 

(7a) 

sin  0ext  =  (tfint  +  47TM0)sin  0int. 

(7b) 

We  note  that  Eq.  (4)  for  an  obliquely  magnetized  nano¬ 
contact  formally  coincides  with  the  equations  derived  for  the 
case  of  normally9  and  in-plane 15  magnetized  nanocontacts. 
The  only  difference  is  in  the  explicit  expressions  for  the 
FMR  frequency  <w0,  nonlinear  frequency  shift  coefficient  N, 
and  the  dispersion  coefficient  D.  Therefore,  we  were  able  to 
easily  generalize  the  results  previously  obtained  in  particular 
cases  to  a  general  case  of  an  obliquely  magnetized  magnetic 
nanocontact. 


N=-2M, 


/  2  .  9 

3(x)q  (  3 cousin  0i 


dMr, 


(O0 


To  describe  the  spin  wave  mode  generation  in  an  ob¬ 
liquely  magnetized  magnetic  nanocontact,  we  will  use  the 
standard  Hamiltonian  formalism,17  successfully  used  before 
for  the  analytical  analysis  of  spin-polarized  pheno¬ 
mena.10,1115  Using  this  approach,  one  can  derive  an  approxi¬ 
mate  equation  for  the  dimensionless  complex  spin  wave  am¬ 
plitude  a{t,  r)  of  the  variable  magnetization  normalized  as 


III.  RESULTS  AND  DISCUSSION 

We  solved  Eq.  (4)  numerically,  assuming  circularly  sym¬ 
metric  profile  and  harmonic  time  dependence  for  the  excited 
spin  wave  mode  [i.e.,  using  ansatz  a=a(r)e~“"r  in  Eq.  (4)], 
and  calculated  threshold  currents  corresponding  to  the  exci¬ 
tation  of  a  linear  propagating  spin  wave  mode  (dashed  line) 
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FIG.  1.  (Color  online)  Main  panel:  Dependence  of  the  threshold 
current  Ilh  on  the  out-of-plane  angle  of  the  internal  bias  magnetic 
field  6mt.  The  dashed  line  shows  the  threshold  current  /|“  of  exci¬ 
tation  of  the  linear  mode,  the  dash-dotted  line  the  threshold  current 
4ul  of  excitation  of  the  nonlinear  bullet.  Dashed  vertical  lines  show 
the  critical  magnetization  angle  6a,  separating  regions  of  excitation 
of  the  bullet  mode  (for  An,  <  0cr)  and  linear  propagating  mode  (for 
Ant  >  8CI),  and  “linear”  magnetization  angle  0lin,  above  which  (0int 
>  <?ijn)  the  self-localized  mode  does  not  exist.  The  parameters  are 
4ttM0  =  8  kG,  Hext=10  kOe,  \ex=5  nm,  L= 5  nm,  q:g  =  0.02,  e=0.3, 
and  Rc= 20  nm.  Inset:  Dependence  of  the  external  out-of-plane 
magnetization  angle  0ext  on  the  internal  angle  Ant-  The  parameters 
are  4ttM0=8  kG  and  Hexl=  10  kOe. 

and  standing  self-localized  nonlinear  spin  wave  bullet  (dash- 
dotted  line)  as  functions  of  the  internal  out-of-plane  magne¬ 
tization  angle  Ant  at  a  fixed  magnitude  of  the  external  bias 
magnetic  field  //ext=10  kOe.  These  two  threshold  curves  for 
two  qualitatively  different  spin  wave  modes  are  presented  in 
Fig.  1.  The  inset  in  this  figure  shows  the  relation  between  the 
internal  An,  and  external  0ml  magnetization  angles  calculated 
from  Eqs.  (7a)  and  (7b)  for  Hext=  10  kOe.  The  solid  line  in 
Fig.  1  gives  the  minimum  threshold  current  at  a  given  inter¬ 
nal  magnetization  angle. 

It  is  clear  from  Fig.  1  that  at  a  certain  intermediate  inter¬ 
nal  magnetization  angle  0CI  the  nature  of  the  excited  spin 
wave  mode  abruptly  changes,  and  the  threshold  curve 
Ah(  Ant)»  although  remaining  continuous,  demonstrates  a  char¬ 
acteristic  kink  at  the  transition  point. 

This  abrupt  transition  from  standing  to  propagating  spin 
wave  mode  taking  place  with  the  increase  of  the  magnetiza¬ 
tion  angle  An,  manifests  itself  more  dramatically  in  the  cal¬ 
culated  dependence  of  the  spin  wave  frequency  generated  at 
the  threshold  on  An,  presented  in  Fig.  2.  In  this  figure,  the 
upper  dashed  line  represents  the  frequency  of  the  quasilinear 
propagating  spin  wave  generated  at  the  threshold  [see  Eq.  (1) 
for  |a|-0],  while  the  lower  dash-dotted  line  represents  the 
threshold  frequency  of  the  nonlinear  standing  spin  wave  bul¬ 
let  [see  the  first  equation  of  Eqs.  (2)].  An  abrupt  jump  in  the 
frequency  generated  at  the  threshold  (of  the  order  of  several 
gigahertz)  takes  place  at  the  point  of  mode  switching  An, 
=  0CI.  In  the  inset  in  Fig.  2,  it  is  shown  that  the  magnitude  of 
this  frequency  jump  increases  with  the  decrease  of  the  nano- 


FIG.  2.  (Color  online)  Main  panel:  Dependence  of  the  frequency 
generated  at  the  excitation  threshold  on  the  out-of-plane  angle  of 
the  internal  bias  magnetic  field  Ant-  The  dashed  line  shows  the 
frequency  of  the  linear  mode,  the  dash-dotted  line  the  frequency  of 
the  nonlinear  bullet.  Vertical  dashed  lines  show  magnetization 
angles  da  and  0[in  (see  caption  of  Fig.  1).  The  parameters  are  the 
same  as  in  Fig.  1.  Inset:  Dependence  of  the  frequency  jump  Au>  at 
An,=  Ar  (see  main  panel)  on  the  magnitude  of  the  external  bias 
magnetic  field  Hext  for  two  different  nanocontact  radii:  Rc  =  20  nm 
(solid  line)  and  Rc= 30  nm  (dashed  line).  All  other  parameters  are 
the  same  as  in  Fig.  1. 

contact  radius  Rc  and,  also,  slowly  increases  with  the  in¬ 
crease  of  the  external  bias  magnetic  field. 

The  calculated  dependence  of  the  critical  internal  magne¬ 
tization  angle  8cr(Hcu)  (at  which  the  mode  switching  takes 
place)  on  the  magnitude  of  the  external  bias  magnetic  field  is 
shown  in  Fig.  3.  It  is  clear  from  Fig.  3  that  the  angle  0CI 
decreases  with  the  increase  of  Hext  and  practically  does  not 
depend  on  the  contact  radius.  It  is  also  clear  that  the  differ¬ 
ence  between  this  angle  and  the  angle  Ain  (upper  dash-dotted 
line  in  Fig.  3)  at  which  the  coefficient  of  nonlinear  frequency 
shift  N  goes  through  zero  and  changes  its  sign  increases  with 
the  increase  of  the  bias  magnetic  field  magnitude. 

So  far,  we  considered  the  behavior  of  a  magnetic  nano¬ 
contact  only  at  the  threshold  of  microwave  generation.  In 
this  case,  changes  in  the  internal  magnetization  angle  0m ,  can 
be  caused  only  by  the  changes  in  the  magnitude  Hext  or  di¬ 
rection  0ext  of  the  external  bias  magnetic  field. 

In  the  above-threshold  case,  when  the  amplitude  of  the 
excited  spin  wave  mode  a  is  substantial,  there  exists  another 
possibility  for  the  variation  of  the  internal  magnetization 
angle  0iM.  The  increase  of  the  mode  amplitude  (or  the  ampli¬ 
tude  of  microwave  precession  of  magnetization)  a  leads  to 
the  reduction  of  the  static  magnetization  in  the  free  layer 
M0=>M,=(l-2\a\2)M0,  and,  due  to  the  boundary  conditions 
[Eqs.  (7a)  and  (7b)],  to  the  changes  in  both  the  magnitude 
//jnt  and  direction  0m ,  of  the  internal  bias  magnetic  field. 

The  decrease  of  the  effective  static  magnetization  M.  in 
the  free  layer  caused  by  the  microwave  precession  will  lead 
to  the  increase  of  the  internal  magnetization  angle  An,- 
Therefore,  if  the  “threshold”  (|fl|  =  0)  operating  point  corre- 
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FIG.  3.  (Color  online)  Dependence  of  the  critical  angle  6CI,  at 
which  transition  from  bullet  to  linear  mode  occurs,  on  the  magni¬ 
tude  of  the  external  magnetic  held  Hext.  Solid  line,  nanocontact 
radius  Rc=  20  nm;  dotted  line,  Rc  =  30  nm.  Dash-dotted  line  shows 
the  angle  0lin,  at  which  the  nonlinear  frequency  shift  vanishes,  N 
— >  0  (this  angle  sets  the  upper  limit  of  existence  of  the  bullet  mode). 
All  other  parameters  are  the  same  as  in  Fig.  1. 

sponds  to  the  region  of  excitation  of  a  bullet  spin  wave 
mode,  0mt<  6CI,  the  increase  of  the  bias  current  and,  conse¬ 
quently,  the  increase  in  the  precession  amplitude  can  cause 
transition  to  the  region  of  excitation  of  a  quasilinear  propa¬ 
gating  spin  wave  mode,  9ml>  0CT. 


As  a  result,  with  the  increase  of  the  amplitude  of  the  bias 
current  /,  one  could  observe  in  experiments  performed  in  the 
nanocontact  geometry  abrupt  jumps  in  the  generated  micro- 
wave  frequency,  similar  to  the  jump  shown  in  Fig.  2. 

IV.  CONCLUSION 

In  conclusion,  we  studied  theoretically  the  nature  of  mi¬ 
crowave  spin  wave  modes  that  could  be  excited  by  spin- 
polarized  direct  current  in  obliquely  magnetized  magnetic 
nano-contacts.  It  was  demonstrated  that  with  the  increase  of 
the  out-of-plane  magnetization  angle,  an  abrupt  switching 
between  two  qualitatively  different  (propagating  and  stand¬ 
ing)  spin  wave  modes  can  take  place.  This  switching  should 
manifest  itself  as  an  abrupt  change  of  the  microwave  fre¬ 
quency  generated  at  the  threshold  of  microwave  excitation. 
We  also  believe  that  the  proposed  mechanism  of  spin  wave 
mode  switching  might  be  responsible  for  the  jumps  in  gen¬ 
erated  frequency  observed  in  some  of  the  recent  experiments 
on  microwave  generation  by  spin-polarized  current.6,7 

ACKNOWLEDGMENTS 

This  work  was  supported  in  part  by  the  MURI  Grant  No. 
W911NF-04-1-0247  from  the  U.S.  Army  Research  Office,  by 
Contract  No.  W56HZV-07-P-L612  from  the  U.S.  Army 
TARDEC,  RDECOM,  by  Grant  No.  ECCS-0653901  from 
the  National  Science  Foundation  of  the  USA,  by  Grant  No. 
3066  from  the  Science  and  Technology  Center  in  Ukraine, 
and  by  Grant  No.  M/226-2004  from  the  MES  of  Ukraine. 


*J.  C.  Slonczewski,  J.  Magn.  Magn.  Mater.  159,  LI  (1996). 

2L.  Berger,  Phys.  Rev.  B  54,  9353  (1996). 

3M.  Tsoi,  A.  G.  M.  Jansen,  J.  Bass,  W.  C.  Chiang,  M.  Seek,  V. 

Tsoi,  and  P.  Wyder,  Phys.  Rev.  Lett.  80,  4281  (1998). 

4W.  H.  Rippard,  M.  R.  Pufall,  and  T.  J.  Silva,  Appl.  Phys.  Lett.  82, 
1260  (2003). 

5W.  H.  Rippard,  M.  R.  Pufall,  S.  Kaka,  S.  E.  Russek,  and  T.  J. 

Silva,  Phys.  Rev.  Lett.  92,  027201  (2004). 

6W.  H.  Rippard,  M.  R.  Pufall,  S.  Kaka,  T.  J.  Silva,  and  S.  E. 

Russek,  Phys.  Rev.  B  70,  100406(R)  (2004). 

7S.  I.  Kiselev,  J.  C.  Sankey,  I.  N.  Krivorotov,  N.  C.  Emley,  R.  J. 
Schoelkopf,  R.  A.  Buhrman,  and  D.  C.  Ralph,  Nature  (London) 
425,  380  (2003). 

8S.  I.  Kiselev,  J.  C.  Sankey,  I.  N.  Krivorotov,  N.  C.  Emley,  M. 
Rinkoski,  C.  Perez,  R.  A.  Buhrman,  and  D.  C.  Ralph,  Phys.  Rev. 
Lett.  93,  036601  (2004). 


9J.  C.  Slonczewski,  J.  Magn.  Magn.  Mater.  195,  L261  (1999). 
10S.  M.  Rezende,  F.  M.  de  Aguiar,  and  A.  Azevedo,  Phys.  Rev.  Lett. 
94,  037202  (2005). 

11  A.  N.  Slavin  and  P.  Kabos,  IEEE  Trans.  Magn.  41,  1264  (2005). 
12 G.  Bertotti,  C.  Serpico,  I.  D.  Mayergoyz,  A.  Magni,  M.  d’ Aquino, 
and  R.  Bonin,  Phys.  Rev.  Lett.  94,  127206  (2005). 

13X.  Zhu  and  J.-G.  Zhu,  J.  Appl.  Phys.  95,  7318  (2004). 

14M.  A.  Hoefer,  M.  J.  Ablowitz,  B.  Ilan,  M.  R.  Pufall,  and  T.  J. 

Silva,  Phys.  Rev.  Lett.  95,  267206  (2005). 

15  A.  N.  Slavin  and  V.  S.  Tiberkevich,  Phys.  Rev.  Lett.  95,  237201 
(2005). 

16Y.  Silberberg,  Opt.  Lett.  15,  1282  (1990). 

I7V.  S.  L'vov,  Wave  Turbulence  Under  Parametric  Excitation 
(Springer- Verlag,  New  York,  1994). 

18  N.  N.  Akhmediev  and  A.  Ankiewicz,  Solitons:  Nonlinear  Pulses 
and  Beams  (Chapman  and  Hall,  London,  1997). 


024437-4 


